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The magnetization reversal in a single molecular magnet (SMM) weakly coupled to an electrode 
with spin-dependent splitting of chemical potentials (spin bias) is theoretically investigated by means 
of the rate equation. A microscopic mechanism for the reversal is demonstrated by the avalanche 
dynamics at the reversal point. The magnetization as a function of the spin bias shows hysteresis 
loops tunable by the gate voltage and varying with temperature. The nondestructive measurement to 
the orientation of giant spin in SMM is presented by measuring the fully polarized electric current in 
the response to a small spin bias. For Mni2ac molecule, its small transverse anisotropy only slightly 
violates the results above. The situation when there is an angle between the easy axis of the SMM 
and the spin quantization direction of the electrode is also studied. 
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I. INTRODUCTION 

Magnetization reversal driven by spin-polarized elec- 
tric current had attracted considerable interests over the 
last decade.— i^i^i^i^ Recent experiments^'^ demonstrated 
that a pure spin current accompanied by no net charge 
current can also be used to reverse the magnetization of 
a ferromagnetic particle. In this setup, the particle is at- 
tached to a nonmagnetic metal wire, in which the chem- 
ical potentials of two spin components are split by using 
the non-local spin injection technique4'^ '^°i^^ By chang- 
ing only the direction of injection current, the sign of the 
splitting can be reversed, leading to the magnetization 
reversal in both orientationsj^i Application of only pure 
spin current is appealing for spintronics devices, because 
it helps to reduce critical currents. Joule heat, noise, etc. 

Meanwhile, another family of intensively studied 
nanoscale magnetic materials, the single molecular 
magnets (SMM)^iiliMii5a6ii7ii8ii9 was reported re- 
cently to be trapped in a typical field effect transis- 
tor geometry, allowing electronic transport measure- 
ment to be performed on an individual SMM with 
great tunability.^"'— )S2,i^ The experiment progresses in- 
spired many transport theories of SMMs, including 
magnetic signatures of SMMs in transport, ^^i^iS^i^ 
Kondo effect ^iSSi^i^iH Berry phase^'^^ full counting 
statistics)^ quantum computing)^ cotunnelingj^ and 
vibrational excitation^i 

Ion spins of magnetic metal in an SMM are inter- 
locked to form a collective giant spin, whose two max- 
imally magnetized ground states orient to opposite di- 
rections due to uniaxial anisotropy, a property that im- 
plies to be a promising candidate for high-density infor- 
mation storage. Therefore, one of the important issues 
is how to manipulate and measure the magnetization of 
an SMM^^I^^^^^ '~ 
and reading qubit encoded by SMM 



port approaches. It has been discussed that spin accu- 
mulation can be induced by charge current)^ Besides, 
it has been proposed that spin-polarized electric current 
injected from ferromagnetic electrodes can be used to 
switch the magnetization of SMM.^i^Siiiii^ 
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FIG. 1: Schematic of energy configurations of our setup to 
manipulate the magnetization of an SMM. The splitting be- 
tween I and I Fermi levels of the only source lead is phe- 
nomenologically denoted as V. Throughout the work the 
middle point of f^l^^ is set at so that nl^^ — ±-j. The 
horizontal lines in the SMM region correspond to resonant 
energies to add an extra electron into the SMM via transi- 
tions from state |0, m) to |1, m ± |) , which is tunable with 
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Gold electrode is used in all the SMM transport ex- 
periments by farf22i2ii2^i2^ which is also among the met- 
als (Au, Ag, Al, Cu) employed in nonlocal spin injec- 
tion technique (for a review, see Ref. [4^ and references 
therein). Therefore, it is natural to expect that the mag- 
netization reversal achieved in ferromagnetic particleSii 
could also be realized in SMM, using the nonlocal spin 
injection. After all, the magnetic moment of SMM is 
much smaller than that of ferromagnetic particle. Al- 
though generating considerable spin-dependent splitting 
of chemical potentials with high efficiency remains a chal- 
lenge, it has been demonstrated experimentally by many 
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other approaches such as the spin Hall effect )^^'^^ the 
spin pumping effect^i^ and incidence of polarized light 
into two-dimensional electron gas.— '^'^'^'^ 

Motivated by the progresses in both the non-local spin 
injection and SMMs, in this work, we demonstrate the- 
oretically that a pure spin current induced by the spin- 
dependent splitting of chemical potential (spin bias) in 
a nonmagnetic electrode^i^i^i^i^ is enough to reverse 
the SMM magnetization, as shown in Fig. [1] where V 
phenomenologically denotes the spin dependent splitting 
of the Fermi levels for t and | electrons in the electrode, 
i.e., ^g^''' = ±V/2. We find that in the context of the 
spin bias: (i) to reverse the SMM magnetization, only 
one nonmagnetic electrode is needed. Neither magnetic 
field'^^ nor magnetic contac t^^'^^i"^^ is required, (ii) Only 
a pure spin current without accompanying a net electric 
current flows in the process of reversal (Fig. [5]), which 
avoids the relaxation of magnetization induced by elec- 
tric current;^ (iii) It sheds a light on the mechanism of 
magnetization reversal from a microscopic point of view, 
and may be extended to mesoscopic magnetic particles 
or films. 

Moreover, we will discuss, in the context of using spin 
bias, several effects not addressed or not clarified in the 
previous literatures on the current-induced magnetic re- 
versal in SMM: (i) By analyzing the transition energy 
spectrum(Fig. [3]), we find that the activation energy at 
which the magnetic reversal starts is determined not only 
by the highestf^i^ but also by the lowest transition en- 
ergy, and is tunable by the gate voltage, (ii) The SMM 
magnetization show magnetic hysteresis loop when scan- 
ning the spin bias back and forth. The hysteresis loop 
can be tuned by the gate voltage and shrinks with in- 
creasing temperature (Fig. [4]). (iii) The avalanche dy- 
namics at the magnetic reversal point is demonstrated 
(Fig. [5]), which supports a microscopic magnetization re- 
versal mechanism, (iv) We show that the ground-state 
orientation of the giant spin in SMM can be read out 
noninvasively, by measuring the charge current through 
SMM driven by a small spin bias (Fig. [7]). (iv) The ef- 
fect of weak transverse anisotropy is considered, (v) The 
situation when there is an angle between the SMM easy 
axis and the spin quantization direction of the electrode 
is discussed. 

The paper is organized as follows. First, we will 
show that the reversal mechanism is irrelevant to specific 
model used in Sec. |TT1 The model and general formalism 
of theoretical approach will be introduced in details in 
Sec. mil In Sec, IIV| we present the numerical simula- 
tions of the hysteresis loops tunable by the gate voltage, 
the avalanche dynamics at the magnetic reversal point, 
and the nondestructive detection to the orientation of 
the giant spin. In Sec. |Vl the correction by the trans- 
verse anisotropy is considered. In Sec. IVIl the case when 
SMM easy axis is not coUinear with the spin quantiza- 
tion direction of the electrode is investigated. Finally, a 
summary is presented to compare the advantages of the 
present work to the existing proposals. 



II. MODEL-IRRELEVANT MAGNETIZATION 
REVERSAL MECHANISM 
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FIG. 2: Left: The two branches of molecular many-body 
states considered in the simulations of this work. The ba- 
sic parameters are given in Sec. IIII Bl Vg = —20 mV. Arrows 
indicate all the steps required to reverse the SMM magnetiza- 
tion from —21/2 to 21/2. Right (four panels): The schematic 
of a single step that increases the SMM magnetization by 1. 
This mechanism is supported by the simulation shown in Figs. 
[5] and El 



By far, many models are proposed to describe SMM, 
with extra electrons added into it. These models include 
the giant spin model )^°'^^i^^'^''i'^^''^^ multi-ion modelj^ 
and those based on the density-functional theory i^^i^^ 

Despite the model employed, one can always select out 
two branches of many-body states of SMM. One is for the 
neutral SMM, the other is for when the SMM is charged 
with an extra electron. Assume the total angular momen- 
tum of the ground state of the neutral SMM is S, which 
has 25-1-1 states for different z component of total angu- 
lar momentum, denoted as |0,to) (m £ [—S, S]). Adding 
the extra electron is like coupling two angular momenta, 
leading to two possible ground-state total angular mo- 
menta for the charged branches, S'i i, respectively. For 
simplicity, we assume the total angular momentum of the 
ground state of the charged branch is 5* -I- ^ , which has 
25" -I- 2 states, denoted as |1,to) {m g [-5 ~ + i]). 
Because of the easy axis anisotropy, the ground states of 
the two branches are |0, ±5) and |1,±(S'-|- i)), respec- 
tively. As an example, two such branches with uniaxial 
anisotropy barriers are shown in Fig. [21 using the giant 
spin model proposed by Timm and Elstc."^^'"^^ 

Suppose one wants to reverse the giant spin originally 
orienting antiparallel with the z axis, i.e., at the state 
1,-21/2) or |0,— 10), to parallel orientation, i.e., the 
state |1, 21/2) or |0, 10). By connecting SMM to the lead, 
one can generate a sequence of transitions that charge 
SMM with spin-up electrons and discharge SMM with 
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spin-down electrons 

|1, -21/2) ^ |0, -10) ^ |1, -19/2) ^ ... ^ |1, -11/2) 
^ ... ^ |0,0) ^ ... ^ |1, 11/2) ^ ... ^ |1, 19/2) 
^10,10)^11,21/2), (1) 

as shown in Fig. [2l During each of these transitions, the 
magnetic moment of an electron spin is transferred to 
SMM, by either adding an spin-l electron into or remov- 
ing a spin- 1 electron from SMM. The right four panels in 
the Fig. [5] depict a typical step, in which an f spin en- 
ters, flips to I spin owing to the exchange coupling while 
increases the giant spin orientation by one unit, and es- 
capes from SMM. Such a step repeats until the giant spin 
orientation is reversed. 

Energetically, to generate this charging-discharging se- 
quence, the Fermi level for the spin-up electrons in the 
lead should be higher than all the transition energies of 
adding an spin- up electron from a neutral state |0, m) to 
a charged state |l,TO-f i), i.e., 

> ^|l,m+i) - E\0.ni), (2) 

while the spin-down Fermi level of the lead should be 
lower than all the transition energies of adding an spin- 
down electron from a neutral state |0, m) to a charged 
state |1, m — i), i.e., 

< -^ll^m-i) - E^Q.m)- (3) 

Because of the anisotropy of SMM, the spectrum of the 
transitions for all the possible m has a finite width, as 
shown in Fig. [31 Therefore, the splitting of /x-f and 
i.e., the spin bias, must be large enough to overcome this 
spectrum width. This thereby defines a threshold spin 
bias for the reversal. 

In the process of the reversal, only "f electrons enter 
SMM while only J, electrons leave SMM at almost the 
identical rate. As a result, a nearly pure spin current, 
instead of electric current, flows between the lead and the 
SMM. Once the reversal is accomplished, i.e., |1, 21/2) is 
occupied, no more t (i) electron can enter (leave), and 
the pure spin current decays to zero. The magnetization 
reversal mechanism discussed above is supported by the 
simulation results shown in Fig. [S] 

Note that the neutral and charged branches employed 
in the above reversal mechanism universally exist not 
only within the giant spin approximation )^°i^^i^^i'^^i^^i"^^ 
but also in the multi-ion model®" and density func- 
tional theoryjS,v62 jj^ multi-ion model of Lehmann 
and LosS)^ the case of ferromagnetic inter- ion interaction 
corresponds to the type of SMM discussed in the present 
work. The green dash-dotted lines in Fig. 1(a) of their 
paper describe the transitions between the neutral and 
charged branches with the difference of the total angular 
momentum by 1/2. Recent density-functional theory cal- 
culation also concludes that the total angular momenta 
for the ground-state neutral and anionic branches are 



10 and 21/2, respectively^S. Therefore the above rever- 
sal mechanism should be universally described by most 
models proposed by far. 

Later we will also show that the above mechanism is 
not qualitatively affected by the transverse anisotropy 
(Sec. |V| and non-coUinearity (Sec. IVI|) . 



III. GENERAL FORMALISM 
A. Model 

In this work, we want to focus on the possibility of us- 
ing spin bias, and what we need are one neutral branch 
and one charged branch as discussed above. Specifically, 
we adopt the model proposed by Timm and Elste^i^ to 
describe SMM, which consists of the lowest unoccupied 
molecular orbital (LUMO) and the phenomenological gi- 
ant spin (GS) S of the molecule, 

'^e{T,i} 

-{D + 5DY,na)Sl + H' + E,,^{Vg) 

a 

-Js ■ S, (4) 

where the first two terms depict the LUMO, ria- = rfj-rfo-, 
and daidl) are the annihilation (creation) operators for 
the LUMO, whose on-site energy is tunable by a gate 
voltage Vg. — e is the electron charge. U is the on-site 
Coulomb repulsion. 

The third and forth terms are for the anisotropy of GS, 
where D describes the easy axis anisotropy and 6D ac- 
counts for the correction to the easy axis anisotropy by 
the occupation of LUMO. H' describes the transverse 
anisotropy. We formally consider its possible leading 
terms 

H' = B2{Sl + S^_) + Bi{SX + St). 

(5) 

For Mni2ac, there is usually only term. We also in- 
clude i?2 term because we want to investigate the effect 
of the transverse terms in a general way. For Mni2ac, B^ 
are several orders smaller than the easy axis anisotropy 
D (Table 4.1 of Ref,19). Therefore, we assume that the 
extra electron brings no correction to them. Following 
most experiments and theories by far ,^°i^"'^'^^i'^'^'^^i®° we 
assume the extra electrons change only the magnitude of 
the easy axis anisotropy. 

The fifth term is the energy of ions that form the giant 
spin, which is also proportional to Vg in the same way as 
LUMO. The last term describes the Hund's rule coupling 
J between the giant spin S and the electron spin s in the 
LUMO. s = ^^d)a{~^)apdi3, where ~& are the vector of 
Pauli matrices. 
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The Hamiltonian for the electrodes used to probe SMM 
reads 

^^Icad = ^ ^kacl^^Ckar , (6) 
k.a,T 

where cl.^^{ckaT) is the creation (annihilation) operator 
for a continuous state in the a{E {S, D}) lead with energy 
Efcc and spin r G In reality, there should be 

an angle 6 between the easy axis of SMM and the spin 
orientation in the electrodes, so the spin quantization 
in SMM is denoted a G {Tii}: and in the lead as r G 
{+,—}. The operator for spin-r electron in the lead is 
related to spin-cr operator through an SU(2) rotation. 

The tunneling between the LUMO and the electrodes 
is described by 

6 

(7) 

where 9 e [0, 7r/2]. We set the easy axis of the SMM 
as z axis. For 9 e [tt/2, n] one just reverses the positive 
direction of z axis. In short, the total Hamiltonian for 
the system we are studying is 

-fftotal = ^SMM + ^^lead + H^- (8) 

We believe that the Hamiltonian employed in this work 
captures the physics required for the magnetization re- 
versal, though it is a simplified model. The physical pic- 
ture of this model can be understood as follows. In a 
Mni2ac molecule, eight spin-2 Mn'^"'' ions and four spin- 
I Mn'*^ ions are exchange-coupled. As a good quantum 
number, their total angular momentum may adopt many 
values, referred as different branches. The branch with 
the lowest energy consists of 21 states with a total angu- 
lar momentum S = 10. The branches with other values 
of total angular momentum are much higher in energy. 
Due to the anisotropy along the easy axis, two degener- 
ate ground states of 5' = 10 branch are the states with 
z component of the total angular momentum ^2 = ±5, 
respectively. In this sense the system is regarded as a 
giant spin of 5 = 10, and is simply described by the term 
—DS^. By adding an extra electron, the energy of the 
molecule changes in several aspects: (i) The first is the 
on-site and charging energy to add this electron. This 
is described by tQ^^Ua- If we set the Fermi level of 
the lead as the reference point, this energy can be com- 
pensated by the gate voltage, so it is absorbed into the 
term (eo — eV^) X^o- "cr- '^^'^ ^ term is due to adding 
the second excess electron to the same LUMO. To sim- 
plify the problem, we assume the second electron will 
also occupy the same LUMO, and exclude the possibihty 
to occupy other states. This is purely theoretical sim- 
plification, and is believed not to affect qualitatively the 
physical consequences in the present problem. We have 



to emphasize that the energy of the ions i?ion(Vg) is also 
tunable to the gate voltage Vg in the same way as the ex- 
cess electron. It does not have to explicitly appear in the 
Hamiltonian because it is counted in the energy of each 
many-body state of the SMM, thus can be discarded, 
(ii) Second, the added electron will interact with the gi- 
ant spin to form spin 5 -f 1/2 and S* — 1/2 branches of 
states. The energy difference between these two branches 
can be characterized by the term — JS • s, where J can be 
found by the splitting between the two branches by us- 
ing first-principle calculation^Si because of the splitting 
is around 2JS. (iii) The third is the correction of the 
anisotropy due to the excess electron. Because the an- 
gular momentum of electron spin is much smaller than 
that of the giant spin, the variation of anisotropy in the 
presence of the excess electron, which mainly leads to a 
curvature change in the energy spectrum, can be roughly 
absorbed into the correction parameter 5D. 

Either the giant spin model proposed by Timm and 
Elste^i^^ or by Romeike and co-workers22i2^ describes 
the above physical picture. Both characterize the 
many-body eigen states by the electron occupation, the 
total angular momentum (note that u = ± corre- 
spond to the total angular momentum — =p i), and 
the z component of total angular momentum. When 
calculating the matrix elements {i\dc\i), the Clebsh- 
Gordan coefficients^ correspond to the linear combi- 
nation coeflicients^l a^,/3^. Most importantly, both 
model are able to capture the main features of the ex- 
periments, e.g., the sophisticated magnetic excitations 
and the negative differential conductance observed in the 
experiments. ^'5 

Besides, Lehmann and Loss^ think that the inclusion 
of the excess electron with respect to the uncharged SMM 
should start with a multi-ion model, in which TV ion sites 
are considered, each with an ion spin s. Nearest ions are 
coupled, by either ferromagnetic or anti-ferromagnetic 
exchange interaction. The excess electron can occupy 
and hop among any of these ion sites. The Hund's rule 
coupling between the excess electron and each ion is lo- 
cal, as well as the anisotropy. When considering anti- 
ferromagnetic inter-ion coupling, the ground state adopts 
a zero total angular momentum and apparently is not the 
case for Mni2ac (but valid for other SMMs, such as Mn4 
dimer— i^). The spatial selection rules they predicted 
mainly occur for the anti-ferromagnetic case, thus will 
not be considered in this work. When considering ferro- 
magnetic inter-ion coupling, the ground state adopts a 
maximal angular momentum Ns. They considered only 
one excess electron. In this case, the electron is free to 
hop and couple to all the ions. As a result, the local 
Hund's rule coupling and anisotropy give rise to global 
giant spin properties. In a word, the single excess elec- 
tron and ferromagnetic inter-ion case of Lehmann and 
Loss's model shares the same spirit of those by Timm 
and Elste and Romeike et al. 

The spatial selection rules are also predicted by using 
the density functional theory.^^ However, because of the 
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lead 100 nm) used in the experimentSS, is much wider 
than the size of the molecule (< 5 nm), we think the 
spatial selection rules, which need precise contacts be- 
tween the lead and the ion sites of the molecule, could 
be smeared in realistic samples. 

Though the coexistence of electron-phonon interaction 
and magnetic excitation is observed the phonon fre- 
quency is beyond the energy scale of the current work. 
For example, the phonon frequency observed by Heersche 
et al^ is about 14 meV, while the magnetic excitations 
observe d^°i^^ or in this work (Fig. and the spin bias 
(Fig. |4]) are of order of meV. Therefore, we do not con- 
sider electron-phonon interaction and its related effect 
in this work, e.g., Franck-Cordon blockade^ or thermal- 
activated effect induced by spin-phonon interactioni^ 



B. SMM states in absence of transverse anisotropy 
and parameters for simulations 

We will use the eigen states of Hsmm when 62 = 64 = 
as unperturbed states ^i^^ The transverse anisotropy will 
be taken into account by perturbation for small B2 and 

When B2 = B4 = 0, Hsmm leads to four branches of 
states for the isolated SMM denoted by \n,m)'^, where 
n(= 0, 1, 2) is the LUMO occupation and m is the quan- 
tum number for {Sz + Sz), the z component of total an- 
gular moment;^ Degeneracy index 1/ only appears when 
n = 1. The four branches are: 

The empty branch 



|0,m) = |0) 



LUMO 



(g) \m) 



GS, 



(9) 



where m € [— S", S"]. 

The two singly-occupied branches 



|2,m) = |n) 



LUMO 



(?) \m) 



GS, 



(11) 



where m G [— S", S"]. 

One can refer to Fig. 2 of Ref. to have a direct 
impression of these four branches. But different from 
Ref. I41!, in this work the higher two branches are far 
above the lower two branches because of large J and U. 

We adopt the parameters based on recent experiments 
and first-principles calculations for Mni2ac {S = 10) as 
D = 0.056 meV, SD = -0.008 meV,^° and J = 3.92 
meVi^liSi Because eg can be compensated by V^, we 
set eo = for convenience. We choose U = 25meV, 
which is comparable to the width of Coulomb diamond 
in experimentsi^i^ For the above parameters (large J 
and U), the two highest branches |2,m) and |l,m)+ are 



neglected in the present work because the branch |1, to)"*" 
is about 2SJ (about several tens of meV) above the 
branch |l,m)~, and the branch \2,m) is even higher. In 
the following numerical simulations, we consider only the 
branches |0,r7i) and |l,m)~. By choosing suitable gate 
voltage Vg, these two branches can be nearly degenerate 
with respect to the Fermi levels of the leadsj^- A typical 
situation of = —20 mV is shown in Fig. ^ 



C. Perturbative correction to SMM states by 
transverse anisotropy 

For Mni2ac, the transverse anisotropics B2 and B4 are 
several orders smaller than the easy axis anisotropy D 
(Table 4.1 of Ref. [11). For B2 < D/S"^ and B4 < D/S'^, 
they can be taken into account by the standard pertur- 
bation calculation. Note that degenerate states such as 
0, ±1) and |0, ±2) are coupled by H' , so one has to per- 
form a degenerate perturbation calculation. We consider 
the first-order correction to the states and the second- 
order to their energies (please refer to Appendix |B] for 
details). 

In the presence of weak B2 and B4, the eigen states can 
only be approximately labeled by the quantum number 
m of (Sz -\- Sz), and becomes a linear combination of all 
the states with the same LUMO occupation;^ i.e.. 



|1,to)^ = a^,\ i)LUMO«i|TO+i)GS+/3™| T)LUMO®|m-i)GS, 

(10) 

where m £ [-S-^,S+^] for = m £ [-S+^,S-^] 
for 1/ = +. 

The doubly-occupied branch 



|0,m)p 



|l,m)p 



^ C°„,,„,|0,m'), 

m' — — S 

rn' — — S—^ 

5_ 1 



(12) 



where p stands for perturbed states by B2 and B4. One 
can expect that C^"^ ~ 1, i.e., \0,m)p are mainly con- 
tributed by |0,m), and |l,m)~ by |l,m)~. 

The projection of magnetization along z axis for the 
perturbed states are obtained, for the branch |0, m)p. 



- E 



(13) 



and for the branch \l,m)p 



E E l^™',™l'™''(14) 



m' = -S 



and rrip = m when 62 = 64 = 0. 

According to Eqs. (|Bip and (|B2p . the vahdity of the 
perturbation requires that 



« 1, 



(15) 
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for arbitrary i,j G {|0, m), |1, m) }. By using Eqs. (|B3|) 
and (|B4[) . one can estimate that H'j^ can be as large as 
628^ and i?4S''* when m ~ 0, and Ei — Ej can be as 
small as D when m ^ 0. Therefore, the perturbation 
only apphes for B2 < D/S'^ and B4 < D/S^, which 
are also reasonable values for realistic Mni2ac molecules 
(Table 4.1 of Ref. [H). 



D. Pauli rate equations 

When connected to the leads, the eigen states of i?sMM 
can transit to each other by exchanging electrons with the 
lead. In the weak-coupling regime and when neutral and 
charged states are nearly degenerate, the sequential tun- 
neling is dominant. The transitions are well described 
by the Pauli rate equations of a reduced density matrix 
spanned by the eigen states of SMM. In this approach. 
Born approximation and Markoff approximation are em- 
ployed, and iJx is treated by perturbation up to the sec- 
ond order.^ (Please refer to Appendix [A] for details.) 
The rate equation can be expressed in a compact form. 



(16) 



where < < 1 are the probability to find the state 
i. In this work, i or j belongs to the 25 + 1 = 21 states 
from the branch \0,m)p and 2(5-1- ^) + l = 22 states from 
the branch |l,m)p , the off-diagonal and diagonal terms 
of the coefficient matrix are given by 



(17) 
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FIG. 3: (a) Resonant energies to add an extra electron into 
SMM via transitions from state |0,m) to |l,m±i)~ for three 

T / 

different Vg. The middle point of fiy is set as energy zero 

point so that /Xg^"^ = ±V/2. (b) Zoom-in when Vg — —18.8 
mV, where the middle point of the entire spectrum is aligned 
with 0, i.e. the middle point of the Fermi levels fJ.y^. (c) 
Zoom-in of the lowest transitions in the middle panel. The 
notation |1, m ± 5)" — [0, m) is short for E^-^ m± i )- — -Bio.m) • 
Note that -E|i__(„± 1 j^- - £|o,-m) = E^j^^^^^i.^- - £|o,m) are 
degenerate in the absence of magnetic field. Basic parameters 
are given in Sec. IIIIBI 



where 



Ki, = T{\m\mco^'iO/2)f{E,-E,-^^t) 
+ sin\e/2)f{E,-E.,-ti-)] 
+ \{j\d^\i)\^[cos\6/2),f{E, - E, + 

+ Bm\e/2)J{E,-E,+^l-)]}. 

(18) 

where the Fermi distribution f{x) = l/[exp(x/fcBr) + 
1] is spin-resolved. The coupling between LUMO and 
the a lead is assumed to be a constant parameter P^ = 
27r^j. |VfcQ,pi5(ti; — efc„) = P for nonmagnetic lead. One 
just replaces t by J, and exchanges -I- and — to obtain 
P-ttj- Because the drain lead is non-magnetic and not 
subjected to the spin bias, one can assume /i^ = /ijj 
for simplicity. For the source, /ig — //g only when 

e = o. 

All the physical quantities can be expressed in terms of 
Pi, such as LUMO occupation J^i^iPi^ the SMM mag- 
netization '^i'nT'pPi, and the a current flowing from the 
a lead to SMM 



where rii and rrip correspond to the n and nip in the state 
i = |n, m)^. 

Notice that since we consider only the branches |0, m) 
and |l,m)~, the terms such as {i\da\k){j\dl^\i) in the gen- 
eral formalism for sequential tunneling*''^ lead to k = j, 
which remove the off-diagonal terms Pi^j from the equa- 
tion of motion for Pi and from the current formula, and 
lead to the Pauli type rate equations of only diagonal 
terms. (Please refer to Appendix \X\ for details.) 

Although the above Pauli rate equations formalism is 
widely employed for the molecules weakly coupled to the 
electrodes,-^'^'^'22;^i2adailL its validity in this work de- 
serves some discussion. In the weakly-coupled regime, 
i.e., the lead-molecule coupling P is so small that be- 
tween two consecutive electron tunnelings (~ 2/P), there 
is enough time for the molecule to relax to the eigen states 
of i?sMM- For example, suppose the molecule is origi- 
nally at the state |0, m — i) = |0)lumo ® — ^)gs- The 
molecule will transit to the state | t)LUMO CS) |m — 5)03 
after a spin-up electron tunnels in. Note that | t)LUMO 
\iTi — ^)gs is not an eigen state of i?sMM, it will relax to 
the eigen state I l,m)~ = a" | l)-LVMo'S)\m+^)GS+l3~ \ T 



.p. 



(19) )lumo ®\m— ^)gs- The time scale of this relaxation is 
~ 1/J, since it is the term — Js-S that couples the states 
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FIG. 4: (a) Magnetic hysteresis loops for diflerent gate voltage 
Vg, when scanning the spin bias V back and forth. T — O.OIK. 
Arrows indicate the scanning direction of V. The scanning 
is assumed to be slow enough to allow the system relax to 
steady state. The triangle corresponds to the magnetic re- 
versal point for the simulation in Figs. [5] and [6] (b) The 
activation energy vs Vg. It consists of two slopes connect- 
ing at Vg = —18.8 mV. The left slope is determined by the 



orientation using the spin bias in Sec. IIV CI 

As discussed in Sec. IIII B[ our numerical simulations 
are based on the Pauli rate equations for the 43 states of 
branches |0,to) and |l,m)~. The parameters are already 
given in Sec. IIII Bl We numerically solve Eq. (fTB|) using 
Runge-Kutta method with the relative error smaller than 
j^Q-io^ The steady solutions are obtained by relaxing 
the equations until for each state, its relative probability 
change with respect to the last iteration step is less than 
10-7. 

We have checked the results for total 84 states. For 
the energy scale of this work, the higher two branches 
|1,to)"'" and \2,m) are hardly occupied. Therefore, the 
abandonment of these two branches is justified. 



A. Gate voltage tunable magnetic hysteresis loops 

When 9 — and B2 = B4 ~ 0, the rate equations 
expressed by Eq. indicate that, besides the selection 
rules |rt — n'l = 1 and |to — to'| — i, to trigger a transition 
0,to) \l,m + a)^ that adds a CT electron onto LUMO, 
two conditions are required: 

(CI) Energetic requirement must be satisfied so that 



\l,m+a) 



transition energy E\i^± 



19/2>- 



0,±10>) 



the right slope is de- 



termined by the transition energy Ei-^ 



|l,±5/2)- 



E\o. 



±2). 



Cir- 



cles correspond to Vg=-17.2, -18.8, and -19.6 mV in the left < m+a)- ~ E\ 



(C2) The state |0,m) should be occupied. 
Similarly, to trigger a transition |1, m -I- a)~ — *■ |0, m) 
that removes a a electron from LUMO, it is required that: 
(C3) Energetic requirement must be satisfied so that 



panel, (c) Temperature-dependent magnetic hysteresis loops 



0,m) ! 



as a function of V for V„ 



-18.i 



mV. Results for experi- 
mental temperatures 3K (Ref. [20h and 0.3K (Ref. [21I ) are 
presented. Other parameters are given in Sec. IIII Bl 



I i)-LVMO®\m+\)GS, and | T)lumo ® I'ti- |)gs- To make 
sure that the relaxation happens before the next electron 
tunneling event, J ^ T must be satisfied. Therefore, in 
this work we use J = 3.92 meV from the first principles 
calculatiori2LSi and F ~ 0.01 meV. 



IV. MAGNETIZATION REVERSAL AND 
DETECTION USING SPIN BIAS 

In this section, we will present the numerical results 
when = and B2 = B4 = 0. The cases when 
-82,-64 <C D and 6^0 will be presented in Sees. IVj 
and IVII respectively. We will show both cases bring no 
qualitative change to the results presented in this sec- 
tion. When 9 = and B2 = B4 = 0, for convenience, the 

I T / I 

notations become /Xg — > /ig and \n,m)p \n,mY . 

We will first introduce the magnetization reversal in- 
duced by spin bias, including hysteresis loops as a func- 
tion of spin bias in Sec. IIV Al and avalanche dynamics 
near the reversal point in Sec. IIVBI Finally, we will dis- 
cuss the nondestructive measurement to the giant spin 



(C4) The state |1, m -I- cr) should be occupied. 
To generate the charging-discharging sequence from 
|l,-21/2) through |l,21/2) shown in Fig. [H the en- 
ergetic requirements (CI) must be satisfied for each of 
them. We list all the 42 transition energies Et^-^ rn±^)- ~ 
-E|o,m) in Fig- 131 Note that for the present model 
and parameters the low and high bounds of the en- 
tire spectrum happen to be -E| 1-1-19/2)- — -Eio.iio) and 
-£'|i,±5/2)- — -£'to,±2), respectively. Note that the highest 
and the lowest transition energies may differ from sample 
to sample, but the following qualitative results are unaf- 
fected. The reversal Sz = — 10 ^ 10 requires V large 
enough to satisfy 



T _ ^ 

Ms - > -E'|l,±19/2)- ~ E\Q,±W) 



- ^ <- F 
Ms - < ^|l,±5/2>- 



E, 



|0,±2), 



(20) 



which thereby defines a threshold voltage or activa- 



tion energy4i Fig. [3] indicates clearly that the ac- 
tivation energy is determined not only by the high- 



est E\ 



l,±19/2)- 



-E'|o,±io)r^'^ but also by the lowest 



i?| 1-1-5/2)- — -E'|o,±2)- More importantly, as shown in Fig. 
[3] the entire spectra of all 42 transition energies can be 

T / I 

shifted with respect to /ig by tuning V^, which means 
that the activation energy is tunable by the gate voltage. 
The activation energy as a function of Vg is shown in Fig. 



8 



llljb). The activation energy can be minimized when the 
center of the entire spectrum is aUgned with 0, which is 



about Vo 



-18.1 



— 18.8 mV for the present parameters. When 
e.g., Vg = —19.6, the activation energy is de- 



|0,±10> 



When Vg > 



-18.8, 



termined by £'|i,±i9/2 
e.g 

i?l 1-1-5/2)- — E\Q ±2)- Using similar analysis, one knows 
that the reversal £'2 = 10^— 10 requires 



Vg = —17.2, the activation energy is determined by 



- ^ > F 
T ^ 

Ms "2 ^|l,±5/2>- - 



- E\o.±io) 

E\0,±2), 



(21) 



As a result, the magnetization when sweeping V back 
and forth must exhibit a hysteresis loop. As shown in 
Fig. mja), the hysteresis loop varies with Vg because the 
activation energy is tunable by the gate voltage. When 
considering the broadening of the Fermi surface at higher 
temperatures, thermal fluctuation will activate the mag- 
netic reversal before V reaches exactly the required acti- 
vation energy. As a consequence, the magnetic hysteresis 
loop shrinks when the temperature increases [Fig. [UJc)]. 

Note that the existence of the activation energy and its 
tunability to the gate voltage only results from that the 
spectrum of all the transition energies has a finite width, 
and the spectrum must has a finite width because of the 
anisotropy. Therefore, although in the reality the highest 
and lowest transition energies may differ from sample to 
sample, the above results are qualitatively unaffected. 




FIG. 5: The time-dependent magnetization and the spin- 
resolved currents during the magnetic reversal. Solid and 
dashed lines represent the cases in the absence and the pres- 
ence of a weak transverse anisotropy, respectively. The ini- 



tial state is Piq . 



-10) 



= 1 and Pi 



i7^|0,-10) 



= 0. r = o.oiK, 



Vg = -18.8 mV, V = 1.8 meV [ the tri angle in Fig. (Ha)]. 
Other parameters are given in Sec. IIII Bl The positive sign of 
current stands for flowing from source to SMM, and negative 
for from SMM to source. 



Nondestructive detection to giant spin 
orientation 



B. Avalanche dynamics at magnetic reversal point 



We investigate the dynamics at the reversal point 
marked by the triangle in Fig. Hlja). On the left 
side of this point, the activation energy is deter- 
mined by 1-1-19/2)- — ±10), and V is smaller than 



2(i?|; 



-19/2)- 



£-10,-10)) while larger than the rest 40 



transition energies. Therefore, the reversal is blocked at 
the state |0, — 10) only because [1,-19/2}^ can not be 
occupied. Once |1, — 19/2)~ is occupied when V exceeds 
2(i?|i _ig/2)- —i?| 0,-10)), an avalanche of the rest 40 tran- 
sitions will be triggered. In Figs. and El we demon- 
strate the numerical simulation of this avalanche by 
showing the time-dependent probabilities for the states 
of the branches |0,to) and |1,to)~ when Vg = —18.8 mV 
and V — 1.8 meV > 2(i?|i _ig/2)- — i?|o,-io))- As we see, 
all the intermediate states between |0, —10) and |1, 21/2) 
show a clear time sequence. Each of them is first occu- 
pied, then reaches a maximum, and finally decays to zero. 
The time scale of the avalanche process can be estimated 



by ~ 300 X f , only - 
fitting parameter— F 



6 ns if we choose the experiment 
= 8GHz w 0.033 meV. 



A scheme to read out the SMM magnetization has been 
proposed by Timm and Elste,'^'' in which however the 
readout is accompanied by the decay of SMM magneti- 
zation. Here we propose a detection scheme, in which 
the giant spin of a SMM has already been prepared to 
be at either S or —S ground state. Our goal is to detect 
at which of the two orientations the giant spin points, 
and most importantly, without destroying the giant spin 
orientation, by means of spin bias.— We have shown in 
Fig. [5] that, in the presence of only one lead only a burst 
of pure spin current can flow only during the magnetic 
reversal, i.e. one lead is not enough to maintain a steady 
current. Thus, we introduce two leads for this detection, 
as shown in Fig. [Tl^b). A small spin bias V = 0.05 meV 

is applied to only the source lead so that fj^J^ = i-j, 

while the Fermi level for the drain lead is set at 
0. By tuning the gate voltage Vg, the Fermi levels of 

T / 1 

leads be located aligned with the two degen- 

erate transition energies i?|i.±2i/2)- ~ -E'|o,±io)- Taking 
Figs. [2ljcl)-(c4), for example, the giant spin orientation 
is initiated at = -10, i.e., either |1, -21/2) or |0, -10) 
is likely to be occupied. On the left side of Figs. [2{c2)- 
(c4), the molecule is at the state |0, — 10); on the right 
side of Figs. [Ii;c2)-(c4), the molecule is at [1,-21/2). 
Notice that the transition between them involves only J, 
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FIG. 6: During the magnetic reversal, the time-dependent 
probabihties of the branches |0,m) and |l,m)~. Thick and 
thin lines represent the cases when B2 = B4 = and 10^52 ~ 
10^54 — D, respectively. Note that the states from jO, —2) 
through |0, 2) are explicitly reshaped by the B2 and B4. All 
the parameters are the same as Fig. [S] 



electron. As shown in Fig. [7]Jb), when 

Ms < -E'|l,±21/2)- - -E^|0,±10) < Md' 



(22) 



a spin-down electron can be injected from the drain lead, 
inducing the transition |0,--10) — > 21/2), then leaks 
to the source lead, making the molecule recover to the 
state |0, — 10). Such process repeats continuously, lead- 
ing to a steady fully polarized spin-down electric current 
flowing from the drain to the source lead. Since we define 
the direction of current from source to drain as positive. 
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Nondestructive detection to the giant spin orien- 



(a) Transition energies E.-^ 



|l,m±^)- 



E\0,: 



when 



Vg = — 18.8mV. (b) Zoom-in of (a) near -B|i,±2i/2) — -E|o,±io>- 

T / 

By tuning the gate voltage Vg, the Fermi levels of leads /^g/^ 
can be located aligned with the two degenerate transition en- 
ergies -E|i__|_2i/2>- ~^'|o,±io>- [(cl) and (dl)]: Schematics of 
how fully polarized electric current flows when Sz ~ ±10. 
[(c2)-(c4) and (d2)-(d4)]: For Sz = ±10, numerical results 
of source-to-drain current, SMM magnetization, and LUMO 
occupation as functions of Vg in the presence of a small 
V — 0.05meV. T=0.1K. Other parameters are given in Sec. 
ElBl 



the spin-| current is negative, as shown by the negative 
peak in Fig. nic2). 

For the same situation shown in Fig. [7{b), if the giant 
spin is oriented along Sz = 10 as shown in Figs. [DJdl)- 
(d4), only the transition between |l,21/2) and |0, 10) is 
possible, which will generate a positive j current from 
source to drain. Because of the degeneracy of 21/2)- 



0,10) 



and En 



|l,-21/2)- 



Ei 



0,-10) J 



the spin polarization 



and flowing direction of the steady current depend only 
on the giant spin orientation. 

Notice that the detection does not change the giant 
spin orientation from initiated Sz — ±10 to other val- 
ues. Specifically, in order to destroy the state initi- 
ated Sz = —10, the state |0, — 10) has to transit to 
|1, —19/2)^ first, which requires two conditions simulta- 
neously: (1) > -E|i^_i9/2)- --B|o,_io). (2) |0,-10> 
is occupied. According to Fig. [31 to make ^ig/Q > 

-E'|i,-i9/2)- — -E'|o,-io)j must be >-18.8 mV. However, 
when Vg > —18.8, only |1,— 21/2) is allowed to be and 
fully occupied, instead of |0, —10), as shown by the occu- 
pation and magnetization in Fig. [7]Jc3) and (c4). There- 
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fore, our detection scheme is nondestructive. 



EFFECT OF TRANSVERSE ANISOTROPY 
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FIG. 8: The energies of the perturbed branches [0, m)p and 
\l,m)p as a function of perturbed projection nip. Circle and 
triangle: B2 = 10"^D, B4 = W'^D. Solid and dashed lines: 
B2 = Bi = 0. Other parameters are the same as Fig. [2] 
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magnetization along z axis rUp [defined by Eqs. ([T3)l and 
([T^ ]. The second is the weak violation to the spin se- 
lection rules^^ according to Eq. so that transitions 
between states jm^ — m'p\ > 1/2 now are possible. 

In Fig. [51 we show the energy as a function of the pro- 
jection of magnetization along z axis mp for each SMM 
states, when B2 = 10^^-D and B/^ = ib^^D. The case 
when i32 = -B4 = is also plotted for comparison. As 
we see, the changes in energy and rUp are ignorably small 
for those states with |?7ip| ~ S. However, the states with 
\mp\ ^ are greatly reshaped, and in some extent mixed 
together. 

The influence of this mixture on the reversal dynam- 
ics is demonstrated by the thin lines in Fig. [51 Let 
us focus on the subfigures from |0, — 2)p through |0,2)p. 
Due to the mixture of |0, -2)p and |0, 2)p, once SMM 
evolves to the state |0, — 2)p, there is certain probability 
that SMM continues evolving to |0,2)p directly, with- 
out through |l,-3/2)p ^ |0,-1)- -> |l,-l/2); ^ 
|0,0)p ^ |l,l/2)p ^ |l,l)p ^ |l,3/2)- ^ |0,2)p. As a 
result, the probabilities of these intermediate states be- 
come smaller compared with those when B2 = = Q. 
Accordingly, the reversal time is also slightly shortened, 
as shown by the dashed lines in Fig. [5l 

The hysteresis loops in the presence of the trans- 
verse anisotropy are shown in Fig. [51 The hystere- 
sis loop when Vg = —19.6 is not changed, because it 



is determined by E^^ ^ 



19/2>f 



0,±10)p 



where both 



l,±19/2>- 



and i?|o^io>p ^re hardly affected by the trans- 
verse anisotropy. On the other side, the hysteresis loop 
17.2 is slightly modified because it is de- 



when Vg — 
termined by E,^ 



|l,±5/2>p 



E\o. 



±2>p 



while E, 



l,±5/2) 



- and 



E\ 



o,2)p are reshaped by the transverse anisotropy. 

According to Fig. [H the energies of |0,±10)p and 
|l,±21/2)~ are barely affected by weak B2 and B^. 
Therefore, the measurements discussed in Sec. II V CI are 
not affected noticeably. Although the spin selection rules 
now allow transitions between |1, 21/2)~ and states other 
that |0,10)p, e.g., |l,21/2)- ^ |0,9)p, these kind of tran- 
sitions have probabilities of order of {B2/D)'^ ^ 10^*, 
thus is hard to be measured. 



FIG. 9: The same as Fig. H except B2 = W'^D and B4 = 

In this section, we will show that the small transverse 
anisotropy only weakly affects the three main results pre- 
sented in Sec. IIVI and brings no qualitative change. 

The transverse anisotropy leads to two main correc- 
tions of the SMM eigen states. The first is the quanti- 
tative correction to the energy and the projection of the 



VI. NON-COLLINEAR CASE 

In this section, we will analyze the influence to the 
three main results presented in Sec. IIVI when considering 
an angle 9 between the easy axis of the SMM and the spin 
quantization direction of the source lead. Because we 
have shown that the weak transverse anisotropy brings 
ignorable effect in Sec. [Vl we only consider B2 = B4 = 
in this section. 
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FIG. 10: (a) In the presence of an angle Q between the easy 
axis of SMM and the spin quantization direction of the lead, 
there are four possible injection and leakage processes, marked 
by the arrows. cos^(&/2) or sin^(6'/2) indicates the relative 
probability of each process, (b) The reversal time will be 
prolonged by l/cos(S). [(c) and (d)] Magnetization and spin- 
resolved currents for the same situation shown in Fig. [5] for 
different Q. 



A. Reversal dynamics 5 7^ 

When 0^0, the + (-) electrons in the source lead 
can be injected into SMM as t(i) or KJ.) electrons, with 
probabilities of cos^ | and sin^ |, respectively. This will 
prolong the reversal time. We will illustrate it by con- 
sidering the same situation of Sec. IIVBI As shown in 
Fig. llOf a). when \ electrons are injected from /ig with 
the tunneling rate oc cos^ |, they are also leaking to /ig 
with the tunneling rate oc sin^ |. From the viewpoint 
of SMM, the t electrons are injected from the lead at 



2 e 



cos( 



an equivalent tunneling rate cx cos — sin 
Similarly, an extra ratio cos is also subjected to the J, 
electron leaking to the lead. As a result, the reversal 
time will oc l/cos6', as shown by Figs. fTUT c) and (d). 
As shown by Fig. fTOTb). Q should be at least > 0.457r to 
increase the reversal time by one order. In other words, 
the reversal time will not be prolonged dramatically un- 
less Q is very close to 7r/2. Of course, the reversal will fail 
when Q = 7r/2, i.e., when the spin quantization direction 
of the lead is perpendicular to the easy axis of SMM. Be- 
sides, as shown by Fig. fTUTd). during the reversal, there 
is still a pure spin current flowing even in the presence 



of 6' 7^ 0, although 9 will introduce the extra injection of 
I electrons from and the extra leakage of \ electrons 
to /.tg . This is because these two effects will cancel with 
each other since both of them are proportional to sin^ | . 



B. Non-destructive detection 7^ 



(a) - 
H — 










T;i 




0.0 
-0.1 
-0.2 



(Cl) 


V 








Sz=-10 




e=7i/4 







(di) 






Sz=10 




e=7i/4 









^0.0 



3 

o 



0.1 
-0.2 

0.1 
0.0 

-0.1 



(c2) 






Sz=-10 
9=371/8 



(d2) 


1 Sz=10 

1 e=3ji/8 


V 





(c3) 


1 Sz=-10 
\ 6=71/2 


\ 





(d3) 


1 Sz=10 
\ e=7:/2 


\ 





0.2 
0.1 
0.0 

0.2 
0.1 
0.0 

0.1 
0.0 
-0.1 

18.4 18.8 19.2 18.4 18.8 19.2 
Vg [mV] 



FIG. 11: (a) Current tends to flow from drain to source when 
-E'|i,±2i/2>- ~ -E'|o,±io) is between /Xg and /ij/^. (b) Current 
tends to flow from source to drain when 2?^ ±21/2)- ~^|o,±io> 
is between /ig and /ij/^. [(cl)-(c3)] The same measurement 
scheme as Fig. ^cl) except for 61 7^ 0. [(dl)-(d3)] The same 
measurement scheme as Fig. [7jd2) except for S 7^ 0. 

As for the non-destructive measurement in Sec. IIV C[ 
the qualitative nature that spin-up (-down) current is 
favored when Sz = S {—S) is irrelevant to 6. Besides, our 
numerical results show that the measurement is still non- 
destructive. However, in the presence oi 9 0, both | 
and I electrons can tunnel via the Fermi level /ig and /ig . 
As a result, the current can flow along both directions, 
different from the 9 = results shown in Figs. [2tc2) and 
(d2). 

We will use the case Sz = —S to illustrate this dif- 
ference, as shown by Fig. [TlTcl)-(c3). When Sz = —S, 
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SMM still favors j current to tunnel through it via the 



transition energy £'|i.±2i/2) 
Fig. fTlT a). when 



o,±io)- 



As shown by 



Ms < -E'|l,±21/2>- 



£-10, ±10) < ^ir/"^ < A^s ' 



(23) 



current will favor tunneling from drain to source with 
a relative probability cx cos^(0/2), and this corresponds 
to the negative current on the left of Figs. [TlTcl)-(c3): 
while as shown by Fig. IllT b). when 



A'S < A'd ^ < -E^|l,±21/2)- - -E|o,±lo) < Ms ; 



(24) 



current will favor tunneling from source to drain with a 
relative probability oc sin^((?/2), and this corresponds to 
the positive current on the right of Figs. [TlT cl)-(c3). 

The case for 5^ = 5 is similar, except that only f 
current is flowing, and the relative ratio between the 
magnitude of the negative and positive currents becomes 
sin^ I / cos^ I . So one can still use this difference to dis- 
tinguish the orientation of the giant spin, unless 9 is very 
close to 7r/2. When 9 = 7r/2, the negative and positive 
currents are the same, as shown by Figs. [TlT c3) and (d3). 
This is reasonable result, because when the spin orienta- 
tion of the source lead is perpendicular to the easy axis 
of SMM, one should have symmetric results for ] and J,. 



C. Summary for 9^0 

We briefly summarize the influences oi 9 ^ Q. Since 
9 does not change the energy of many-body states, the 
steady-state solutions of hysteresis loops shown in Fig. 
|4] are not affected except 9 ^ ti/2. It will prolong the 
reversal time by l/cos6', as shown in Fig. [TOl There 
are still a pure spin current flowing during the reversal 
when 9^0. The non-destructive measurement shown 
in Fig. [7| is still non-destructive when 9 0. However, 
the measurement signals will be reshaped, as shown in 
Fig. [m Fortunately, one can still employ the relative 
ratio between current flowing along opposite directions 
to distinguish the orientation of the giant spin unless 9 ~ 
tt/2. 



(a) 



two normal leads 



(b) 



two anti-parallel ferromagnetic leads 



(c) 



4. 



two parallel ferromagnetic leads 



(d) 



-Pi 



one normal lead with spin bias 



FIG. 12: [(a)-(c)] Energy schemes employed by the previous 
authors and (d) in this work, (a) Nonmagnetic leads (Ref. 
[STI ). (b) ferromagnetic leads with anti-parallel polarizations, 
and (c) ferromagnetic leads with parallel polarizations (Ref. 
fill ): (d) nonmagnetic lead with spin-dependent splitting of 
chemical potentials. 



reversaL^^' Because no excess spin is transmitted, two 
parallel aligned ferromagnetic leads^ are equivalent to 
two nonmagnetic leads when the charge bias is large 
enough to cover all the transition energies, as shown in 
Fig. [mjc). Moreover, the electric current induced by 
the charge bias between two nonmagnetic electrodes may 
lead to the decay of SMM magnetization.'^^ In contrast, 
the current induced relaxation and low efficiency can be 
avoided in our one-lead model, as only a pure spin current 
flows during the reversal and no lead-SMM electron ex- 
change is permitted when there is no reversal occurring. 
Technically, the ferromagnetic lead can already be at- 
tached to the single molecule with the charge bias voltage 
easily exceeding the required threshold voltage, while 
large spin bias over 1 meV still remains an experimental 
challenge. 

In conclusion, the spin bias or spin current can be 
applied to control and measure the magnetization of 
SMM efficiently 



VII. SUMMARY AND DISCUSSION 

Before ending this paper, we compare the advantages 
of the previous works on current-induced reversal of SMM 
magnetization to our proposal. For the setup with one 
nonmagnetic and one ferromagnetic lead'^^ or the setup 
with two magnetic leads that are not fully polarized)^ 
electrons injected from one lead can leak to the other 
lead, while electrons hopping to one lead can be refilled 
by electrons from the other lead, as shown in Figs. [UKa) 
and (b). The leakage and refilling reduces the efficiency 
of magnetization reversal because only the excess trans- 
mitted spins contribute to the process of the magnetic 
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APPENDIX A: THE RATE EQUATIONS 

In this appendix, we will present the deduction of the 
rate equations Eq. (|16p following the approach intro- 
duced in Ref. ,67 1 . Besides, we will explain why only di- 
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agonal terms of the reduced density matrix are employed 
for the present problem. For simplicity, only 9 — case 
will be addressed. The case for 6* 7^ is a straightforward 
generalization. 

Suppose we have found the eigen states \n) of Hsmm- 
By using the completeness \n) {n\ — unity, the Hamil- 
tonian Eqs. ^ and ([7|) can be rewritten as. 



Hsmm ^^En\n){r 



(Al) 



and 



Ht=J2 E(^fc-("l^-l"^)4aJ")("^l) + /i-c- (A2) 



k,a.cr n,m 



The Liouville equation of the density matrix p of the 
entire system is given by 



dtp = -i[H,p], 



(A3) 



where H = -ffsMM + -fficad + Ht- For arbitrary oper- 
ator O, one introduces the interaction picture Oi{t) = 
e^fiotQe-iHot^ where Hq = i?SMM + i^iead, then Eq. ^ 
becomes 

dtpi{t) = -i[HTi{t),pi{t)]. (A4) 
Integrate and iterate Eq. (|A4p for one time 
dtpiit) « -i[HTi{t),pi{0)] 

dt'[HTi{t),[HTiit'),piit')]]. (A5) 



The reduced density matrix pu (t) of SMM is obtained 
from p{t) by taking the trace over all variables of the 
leads, i.e., in the interaction picture, 

Pmi = TrL/O/. (A6) 

Trace out the lead part in Eq. (jASp 

dtPMiit) = -zTrL[-ffTi(t),/5/(0)] 

dt' Til [Hti (t) , [Hti {t'), pi{t')]]. ( A7) 



The above equation assumed that Ht is switched on at 
t = 0. Prior to this, SMM and the leads are uncorrelated 
and the total density matrix is given by their direct prod- 
ucts, 



p{Q) = pm(0)®pl(0) =/5/(0). 



(A8) 



At this point, we follow Fano to make two key 
assumptionsi^ The first is the Born approximation, 
which assumes that the leads have so many degrees of 
freedom that the effects of interaction with SMM dissi- 
pate away quickly and will not react back to any sig- 
nificant extent, so that the leads remain described by a 



thermal equilibrium distribution at constant temperature 
Pl(0) at all time. 



Piit) ^ PMi{t)(E>PL{0). (Born) 



(A9) 



The second is the Markoff approximation, which assumes 
that due to the rapid relaxation in the leads, correlation 
functions of lead electrons decay on a time scale much 
shorter than the SMM dynamics, i.e., the leads do not 
have memory. This allows replacing the correlation func- 
tions by delta-functions in the rates, which are convolu- 
tions of such correlation functions with the reduced den- 
sity matrix in the sequential-tunneling approximation. In 
this context, dtPMi{t) depends only on its present value 
PMi(i), 

PMi(t') PMiit). (Markoff) (AlO) 



Put the two approximations (jA9p and (|A10[) into Eq. 



dtPMiit) = -zTrL[ffTi(i),/5Mi(0)pL(0)] 

dt'Tlx[HTi{t), [HTi{t'),pMlit) PL 

(All) 

What follows is straightforward calculation by putting 
Hti into Eq. (|A11[) . The first term on the right side of 
Eq. (|A11[) vanishes because 

TrL[cL./(i)/5L(0)] = Tr^[ckc.i{t)pLm - 0. (A12) 
Besides, the cyclic property of the trace is used 

TrL[cL./(i)cfc'a'.'/(t')/5L(0)] = TrL[cL,,(t - t')cfc'a'.'PL(O)]. 

(A13) 

Finally, one returns to the Schrodinger picture, and ar- 
rives at the equation for the arbitrary terms Py = 

(*iPA/b), 

i'j' acT 



-lY.J2^^o.{W){j'\dUmiEf-E,+p^j]p,,, 

i'j' Q(7 

i'j' aa 

- ^ E E 10 (/ 14 Ij) [/(s.' - + 

i'j' Ota 

4EEra(/M^U)(^l4N')/(i?.-i?.'-<m'.' 



i'j' QL(7 
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-lY.T.^''c.{J'K\^'){^'\dl\J)fiE,-E,,-f,'^JP,,, 

i'j' aa 



(A14) 



The above equation is exactly Eq. (2) of Ref. [60 by 
replacing i ^ a, j ^ P, i' a', and j' — > 

In this work, we have calculated three physical quan- 
tities, the total magnetization, the LUMO occupation, 
and the current through SMM. The general form of a 
current flowing from lead to SMM obtained by Lehmann 
and Loss can be rewritten using our notation as^*' 

= eReJ2K{f{E^-E,,+^^l){J\4\^'){^'K\^) 

ii'j 

-f{E,-E,-^,'^J{J\d,\^'){^'\dl\^)}P,,. 

(A15) 



As discussed in Sec. IIIIBi we neglect two high-energy 
branches |1,to)"'" and \2,m), and consider only two low- 
energy branches |0,77i) and |l,m)^. In this context, the 
terms such as {j\dcr\i'){i'\dl^\i) must require that i = j, 
e.g., if \i') = |l,m)~, only (0,m — a\da\l,rn)~ and 
~(1, m|c?J.|0, m — cr) are nonzero, so both \i) and \j) can 
only be |0,m — cr). As a result, Eq. (jAlSP reduces to 

= eJ2K{f(E^-E,,+^^l)\{^y,\^)\' 

ii' 

-i?,-MS)|(z|d,|*')|'}P„(A16) 
which can be further simplified as^ 

II - -e^K -n,)i?f4P,. (A17) 



The LUMO occupation N and the total magnetization 
M involve only the diagonal terms of density matrix. 



N 



cr ijk cr 



(A18) 



and 



M = TrM[/5M(S'2 -I- 



ijk 



ki 



because {sz + Sz) and n^r are eigen operators of 

HsMM- 

Therefore, one only needs to know the diagonal terms 
of the density matrix. By letting j = i in Eq. (|A14p . and 
employing the property {i\da\j'){j'\d'l\i') ^ i — i', one 
can readily show that the equation of Pi only couples to 
other diagonal terms of density matrix. 

Note that the property {i\dcr\j'){j'\dl\i') ^ i = i' 
will be weakly violated in the presence of the trans- 
verse anisotropy B2 and -B4, so that terms such as 
~ (1, m|c?J. |0, TO — (t')p may be also nonzero for \a'\ > \a\. 
However, the relative probability of these terms^^ ^ 
{B2/Df, {Bi/Df, thus can be omitted for ^2,^4 < D. 



APPENDIX B: PERTURBATION OF 
TRANSVERSE ANISOTROPY 

If we denote H^mm = H'^^^+H', the miperturbed eigen 
states of -ff^^-* by "0^"', the first-order correction to tpf'"' 
is given by 



(1) 



E ^j» , (0) 

77(0) _ £;(0) ^3 ' 



E 



(Bl) 



where E^^^ is the unperturbed energy of ipf^' , and iJ'^ — 



(0) 



{ipj^''\H'\ipl'''')- The second-order (the first order is zero 
for the present problem) correction to the energy is given 
by 



/,(0) 



E, 



(2) _ 



E 



E. 



(0) 



E 



(0) 



(B2) 



Finally, the perturbed states and their energies are ob- 
tained as V = ■0^°^ + 0^^^ and E = + E^'^\ Because 
|0, to) and |1, to)^ are not coupled by H' , we will discuss 
them separately. 



1. Branch |0, m) 

The unperturbed energy is given by 



^ = -Dm^ 

|0,m) ' 



(B3) 



(A19) and the perturbation 



(0,m'|i/'|0,m) 

= -52^(5 - m){S + TO + \){S - TO - \){S + m + 2) 6,n'=m+2 
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-B2V{S + m){S -m+ 1){S + m - 1){S - m + 2) (5„/=„^2 
-B4 - m){S + m+ 1){S - m - 1){S + m + 2) 

X y^iS - TO - 2){S + m + 3){S - m - 3){S + m + 4) (S™'=„+4 
-54^(5 + m){S - TO + 1)(S' + TO - 1)(S' - TO + 2) 

X ^(5" + TO - 2)(S' - TO + 3)(S' + TO - 3)(S' - TO + 4) (5™'=m-4 

(B4) 



However, because degenerate states |0,±1} and |0, ±2) 
are coupled by H', Eqs. (|Bip and (|B2[) can not be applied 
directly. One has to perform a linear transformation first 
^-f"^ , = PV'm^ V, so that 

|0,m) ^|0,m)' 



Then one replaces ^ H^°\ H' 



*f ^ in Eqs. 



h' , and V'i°^ ^ 
Bip and (jB2|l to perform the perturbation. 



and 



(*[o;-i)i^'i*[o,\))- 



(*[o?i)i^'i*s?-i)) = 



\ |0,-2) 



l^'l*S)) = (*(o?2)l^'l*;o%) = 0-(B6) 



Under the same linear transformation, iJ'^*'^ and if' be- 



(B7) 



2. Branch |1, m)- 



Because |1,to)^ and |1,to)+ together constitute the 
complete set of the subspace of n = 1, the perturba- 
tion of |1,to)~ has to take |1,to)+ into account. The 
unperturbed states and their energy are obtained by di- 
agonalizing 



/^LUMo(TO+l/2|Gsa i)LUMo|TO + l/2)GS LUMo(m + l/2|Gsa 1^^^°^! T)lumo |m - 1/2)gs \ 

VLUMo(m-l/2|Gs(T i)LUMol™ + l/2)GS lumo(™-1/2|gs(T T)LUMo|m-l/2}Gsy 

/ eo - eVg + 4to + ^ - (P> + <5P>)(to + \f -^J^S{S + 1) - (to + ^){m - i) 

y -y^S{S + 1) - (to + i)(TO - i) eo^eVg-jm+{-{D + SD){m-^)^ 



for TO G {—S* + i, — i}, so that in each subspace of TO, and a|i,s+i/2) = P\i.-s-i/2) = and a|i,-s-i/2) = 

P\l,S+l/2) = 1- 

|1,to)± = «±||)|^+i)+/3±|t)|™-i}, (B9) 



The perturbation is then given by 

'''(l,TO'|ii'|l,TO)' 

= -S.^>^m+2B2[iai,ra-^^{S - to - + to + 1^){S -m-l^){S + m + 2^) 
+ {P^'rP::^\|{S ^m+^){S + m +l)iS-m- ^){S + TO + li)] 

-(5m'=m-2P2[(am')*amy (-S* + W-f i)(5 - TO-j- i)(S' -|- TO - ^){S - TO-|- 1^) 
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HPi'rP''m\j{S + m--){S-m + 1-){S + m -l-){S-m + 2-)] 
-5m'=m+4-B4[«„,)*a^y (5 - m - ^ ) (S* + m + li)(5 - to - li)(S' + m + 2^) 
x^(5-™-2i)(5 + m + 3i)(5-TO-3i)(^ + m + 4i) 
+ - m + + ^ + - m - + to + 1^) 

x\j{S-m-l\){S + m + 2\){S-m-2\){S + m + ?}-)] 

-5,n'=rn-ABi[{atr<n\j {S + m + - TO + + „j _ _ ^ + li) 

xy(5 + TO-li)(5-m + 2i)(5 + TO-2i)(5-m + 3i) 

+(/3rn'')*/3;;y(5 + to - - m + ii)(5 + _ _ + 2!) 

X W(5 + TO - 2]-){S - m + 3i)(S' + to - 3i)(S' - m + 4^)]. (BlO) 



Notice that here to are half- integers, so there is not de- directly, 
generacy problem. One can employ Eqs. (|Bip and (IB2p 
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